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We study quantum chaos in open dynamical systems and show that it is characterized by quantum 
fractal eigenstates located on the underlying classical strange repeller. The states with longest life 
times typically reveal a scars structure on the classical fractal set. 
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At present, the structure of phase space for generic 
classical Hamiltonian systems is well understood both on 
qualitative and quantitative levels. With the increase 
of perturbation's parameter the invariant KAM curves 
are destroyed by the perturbation and are replaced by 
chaotic regions. These regions grow and become inter- 
connected over the whole phase space. In this situation 
the phase space exhibits a hierarchical structure of mixed, 
chaotic and integrable components which continues on 
smaller and smaller scales of the phase space. At larger 
perturbations the measure of stability islands decreases 
until they become negligible for sufficiently strong per- 
turbations. This is a general scenario for the emergence 
of chaos in classical Hamiltonian systems [Q. The well 
known model in which such scenario has been studied in 
detail is the Chirikov standard map 

The investigations of the corresponding quantum sys- 
tems show that the structure of eigenstates is closely 
related to the properties of classical phase space. In 
the integrable regime, eigenstates are located on the in- 
variant KAM curves while in the chaotic regime they 
spread over the whole chaotic component in agreement 
with the Shnirelman theorem |^ . The structure of eigen- 
states can be seen in a pictorial way with the help of 
Wigner function and Husimi distribution [Q . This repre- 
sentation allows to see graphically the qualitative change 
of eigenstates during the transition from integrability to 
chaos. In the chaotic regime they also allow to see scarred 
eigenf unctions in which the probability is concentrated 
near short unstable periodic orbits |^,^. However in the 
chaotic regime, the majority of eigenfunctions are ergodic 
on the energy surface and the energy level statistics is well 
described by the Random Matrix Theory 0| . Here we as- 
sume that the value of the Planck constant is sufficiently 
small and the effect of dynamical localization does not 
alter the properties of eigenstates |^ . 

As a result, in absence of localization, the structure 
of quantum eigenstates in the regime of Hamiltonian 
chaos is now well understood and many concrete exam- 
ples have been studied numerically |^|^ and experimen- 
tally |9j. Conversely, the chaotic objects with fractal 



structure, which appear in dissipative classical dynamics 
such as strange attractors and repellers, were not studied 
in quantum mechanics. The main problem is that in the 
quantum case the dissipation is always accompanied by 
noise and generally one should study the density matrix 
In this way the problem becomes much more 
complicated than the Hamiltonian case and the fractal 
structure of strange attractors/repellers had never been 
seen in quantum eigenstates. 

In this paper we study an open chaotic system in which 
absorption leads to the appearance of a fractal set in the 
classical phase space (strange repeller) . The quantum 
dynamics of the model is governed by a nonunitary evo- 
lution operator. Such unitary breaking, due to absorp- 
tion, was widely used in nuclear and mesoscopic physics 
to describe coupling with continuum and massive leads 
via open channels p^ ]. 

The underlying classical fractal set should affect the 
quantum dynamics and find its manifestations in the 
structure of eigenstates. Indeed it is natural to expect 
that long living eigenstates will be associated with the 
above strange set on which classical orbits live forever. 
Therefore the eigenstates associated with this set should 
strongly influence the scattering process, relaxation and 
ionization into continuum. 

To investigate the quantum fractal eigenstates we 
choose the kicked rotator model with absorbing boundary 
conditions which was introduced in ||l3|| . The quantum 
dynamics is described by the evolution matrix: 



(1) 



where h — —id/ 89, h — 1 and the operator P 
projects the wave function to the states in the inter- 
val [—N/2, N/2]. The quasiclassical limit corresponds to 
A: — !■ cx), r — > with the chaos parameter K = kT = 
const. 

The classical dynamics is described by the Chirikov 
standard map [^j: 



n — n + k sin 



Tn 



+ |(n + rl). (2) 



1 



In this model all trajectories (and quantum probabilities) 
leaving the interval [— N/2] are absorbed and never 
return back. For the classical map the ionization time 
required for a trajectory to reach diffusively the absorb- 
ing boundary is tc ~ = iV^/Z? where I? « is the 
diffusion rate for i^T 1. The independence of tc on N 
requires N/k — const. 

In our study we fixed N/k = 4 and K = 7 so that 
the phase space is completely chaotic with no visible is- 
lands. As a result the classical probability P{t) to stay 
inside the sample [— A''/2, A^/2] decays exponentially for 
t > to- Pit) ~ exp(-7ci) with 7^ = l/t^ = 0.10188. 
The quantum probability follows closely the classical one 
up to the quantum relaxation time scale tq ~ -^/tc/A p^ ] 
where A = 1/A^ is the levels spacing. 

In the classical case, the orbits which are never ionized 
and which stay forever inside the sample form a fractal set 
(strange repeller). To obtain this set we study the evolu- 
tion of M = 2.2 • 10^ classical orbits up to t = 100 map 
iterations. Initially the orbits are homogeneously dis- 
tributed in the phase plane. An example of a set of points 
which are never ionized and form a fractal is shown in 
Fig.l. The fractal nature is demonstrated by magnifica- 
tion of a small part of the phase space which clearly shows 
a structure typical of strange attractors/repellers |p],|l5|. 
Similarly to the problem of diffusion in the Lorentz gas 
[^6|, the information dimension di of this set can be ex- 
pressed via the Lyapunov exponent A w lni^/2 and the 
probability decay rate 7c ~ D/N"^: di = 2 — 7c/ A. In 
our case with 7c « 0.1, A « 1.25 this gives di w 1.92. 

In the quantum system, due to absorption at the 
boundary, all eigenvalues of U move inside the uni- 
tary circle and can be written as A = exp (— ie) = 
exp {—iE — T/2). The imaginary part of e determines 
the decay rate F of an eigenstate. Due to the symmetry 
of the t/-matrix the eigenstates are symmetric or anti- 
symmetric in n and in order to study the statistics of 
F's we restricted ourselves to investigation of symmetric 
states. The general structure of the distribution of F's 
had been found in [|l3|. This distribution dW/dV has 
a gap for small values of F (F < 7c) while for F > 7c 
it drops according to a power law, dW/dr - F-5 for 
tc >> 1 [|l^. A typical example for N ~ 12001 is shown 
in Fig. 2 and confirms the above global structure of the 
distribution. With the increase of the matrix size A^, 
the minimal value of F = Fmm converges to the classical 
value 7c as it is shown in the inset of Fig. 2. The fit of 
numerical data gives AF = 7c — Fmm ~ 6.27c(?~" with 
a = 0.507 ± 0.037. Theoretically we expect that, due to 
fluctuations, AF should be of the order of the distance 
i5 between the eigenvalues of I] in the complex plane. 
Since most of the N eigenvalues e are distributed inside 
the ring of width 7c, then 6 - (7c/A^)^/^ 1/t,. This 
gives AF ~ (5 ~ 1/fg. Therefore, after the fitting of the 
data in Fig. 2, we have 



AF = 7c - Vram ~ 6.27c/ ~/g. (3) 

where g = 7c/ A is the effective conductance of the sys- 
tem. 

To analyze the structure of eigenfunctions in the phase 
space we used the Husimi function obtained from the 
Wigner function smoothed in the intervals An and A0 
(AnAe* = 1/2) (see g). The ratio s = IS.Q / l\n was fixed 
in a way to have optimal resolution. In Fig.l it was equal 
to s = 0.0015(a); 6.25 x 10-4(c);2.25 x 10-4(e) and in 
Fig.3 it was s = 0.0015(a); 0.16(6, c). The Husimi func- 
tions were constructed from antisymmetric eigenstates 
since they had the minimal values of F. The comparison 
between the density distributions in the phase space for 
the classical and the quantum case can be seen in Fig.l. 
The quantum Husimi function for the eigenstates with 
F « 7c reproduces very well the fractal structure of the 
classical strange repeller on very small scales. It shows 
close agreement between the classical and quantum data 
for large N corresponding to a small effective Planck con- 
stant ?ie//- Of course, on very small scales comparable 
with the minimal quantum cell, the quantum density be- 
comes smooth. However, on scales larger than this cell, 
the fractal structure is obvious: we call such states quan- 
tum fractal eigenstates. 

It is natural to expect some analogous of Shnirelman 
theorem for these quantum fractal eigenstates, so that 
in the quasiclassical limit they will be distributed over the 
classical fractal set according to the classical measure. 
The case of Fig.l for the states with F w 7c confirm such 
expectations. However, the situation for the states with 
F = Tmin is rather different (Fig.3). Indeed, generally, 
we observe there the appearance of scars on the underly- 
ing classical strange repeller, the silhouette of which, in 
spite of scars, is still clearly seen. The strength of scars 
grows with decreasing N (increasing heff). However, 
even for the largest N = 59049, the probability distribu- 
tions f(n) projected on the n-axis demonstrate an evi- 
dent difference between the classical and quantum cases 
(Fig. 4). We qualitatively understand this phenomenon in 
the following way: the quantum state with the minimal 
value of F should stay away as far as possible from the ab- 
sorbing boundaries so that quantum interference should 
redistribute probability on the classical set and will lead 
to some scarring around classical trajectories which cor- 
respond to short unstable periodic orbits located near the 
center. Indeed Fig. 4 shows a pronounced peak near pe- 
riod two orbits. There is also some correlation between 
the density of long periodic orbits (Fig. 3d) and quantum 
distributions with not very large N (Fig. 3b, c). However, 
at present, we cannot propose any quantitative explana- 
tions of the scarred eigenstate structure. 

The complex eigenvalues of the evolution operator U 
can be considered as some poles of the scattering matrix. 
Therefore, we can expect that similar quantum fractal 
eigenstates will appear in the problems of chaotic scat- 
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tering in the quasiclassical regime. One of the possi- 
ble models to study such effects is the three discs prob- 
lem where the gap in the F-rates has been discussed in 
pTf . Such quantum fractal eigenstates can be also stud- 
ied in experiments with chaotic light in micrometer-size 
droplets where the classical dynamics is governed by 
a map analogous to the map (2). In conclusion, we have 
demonstrated that the quantum eigenstates can form a 
quantum strange repeller. We conjecture that quantum 
strange attractors, once identified, should have a similar 
structure. 

We thank B.Georgeot for the discussions about 
the Shnirelman theorem for nonunitary operators and 
R.Ketzmerick for useful advises in Lanczos approach to 
the kicked rotator. 
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FIG. 1. Left column: Husimi function of the quantum 
fractal eigenstate for = 59049 and F = 0.1080 7c. 
Right column: classical strange repeller obtained, from ini- 
tially homogeneously distributed orbits, after 100 iterations 
of map (2). The color varies proportionally to the density: 
from blue for zero density to bright red for maximal den- 
sity (the color scale is the same for classical and quantum 
cases on each magnification level). The size of the phase 
region {e,n) is: < 61 < 2%, ~N/2 < n < N/2 for (a) 
and (h); < e < 2tt/3, -N/6 < n < N/6 for (c) and (d); 
0<9< 2tt/9, -N/18 <n< N/18 for (e) and (f). 



FIG. 2. Distribution dW/dT as a function of F for 
TV = 12001. The arrow marks the value of 7c. The inset 
shows AF = 7c — Fmin vs. N in log-log scale. The straight 
line shows the dependence given by (3). To guide the eye, the 
numerical data (points) are connected by lines. 



FIG. 3. Husimi function for quantum eigenstates with 
minimal F = F„i„ for: (a) N = 59049; (b) = 729 and (c) 
A'' = 243. The case (d) shows the classical initial density of 
orbits which survive after 60 iterations forward and back in 
time. The maximum of density (bright red) from (a) to (d) 
is fixed by case (a). 



FIG. 4. Probability distribution /(n) over the unperturbed 
basis n for quantum eigenstates of Fig. la (dotted curve). 
Fig. 3a (dashed curve). The full curve is the classical dis- 
tribution of Fig. lb. Crosses mark periodic orbits with period 
one at (0,0) and period two at (n — ±7r/T for 6 = it/2, 37r/2). 
The asymmetry of /(n) in n for the classical case is due to 
round off computer errors. 
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